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ABSTRACT 

The  plane  elastic  problem  of  a rectangular  orthotropic  region  is 
considered;  subject  to  the  boundary  conditions  of  prescribed  equal  and 
opposite  tangential  displacements  and  zero  normal  displacements  on  the 
upper  and  lower  edges  and  zero  stresses  on  the  remaining  edges.  The 
effect  of  the  stress-free  edges  on  the  stiffness  coefficient  relating  the 
tangential  displacement  and  the  corresponding  shearing  force  is  estimated 
in  the  form  of  upper  and  lower  bounds  for  this  coefficient. 
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Effect  of  Stress-free  Edges  in  Plant;  Shear 
of  a Rectangular  Orthotropic  Region 

INTRODUCTION 

In  what  follows  we  consider  the  plane  problem  of  a rectangular 
elastic  region  which  represents  an  orthotropic  lamina  of  unit  thickness 
in  plane  stress  or  the  cross  section  of  an  infinitely  long  orthotropic  pad 
in  plane  strain.  The  boundary  conditions  consist  of  prescribed  equal 
and  opposite  tangential  and  zero  normal  displacements  on  the  upper 
and  lower  edges  and  zero  tractions  on  the  remaining  two  edges.  Our 
objective  here  is  to  obtain  upper  and  lower  bounds  for  the  stiffness  coef- 
ficient relating  the  prescribed  tangential  displacement  to  the  shear  force 
required  to  produce  this  displacement.  An  elementary  approximate  solu- 
tion to  this  problem  may  be  obtained  by  considering  a state  of  pure  shear 
within  the  region,  which  would  necessarily  exclude  the  effect  of  the  stress- 
free  edges.  In  previous  work  rl],  Read  obtained  bounds  for  the  stiffness 
coefficient  including  the  effect  of  the  stress-free  edges  in  the  case  of 
isotropic  materials  through  the  use  of  the  Prager-Synge  hypercircle 
method  [2].  Our  approach  here  consists  in  the  simultaneous  application 
of  the  principle  of  minimum  potential  energy  and  the  principle  of  minimum 
complementary  energy  to  bound  the  stiffness  coefficient.  These  energy 
principles  have  been  used  previously  to  obtain  bounds  for  influence  coef- 
ficients associated  with  certain  two  dimensional  boundary  value  problems 
in  the  case  of  cantilever  beams  T3]  and  in  the  case  of  cylindrical  shells 
[4].  Our  approach  differs  from  that  in  [ 1 ] in  the  interpretation  of  the 
bounding  functionals  and  generalizes  the  work  in  [ I ] by  considering 
orthotropy  and  by  removing  a restriction  used  in  [1],  that  the 

thickness  to  length  ratio  of  the  rectangular  region  is  small  compared 
to  unity.  The  general  bound  results  obtained  for  the  orthotropic  case 
are  specialized  to  the  isotropic  case  so  as  to  obtain  improved  numerical 
results  for  the  case  considered  in  [ 1 ]. 
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J We  complement  our  bound  calculations  by  presenting  exact 

expressions  for  the  stiffness  coefficient,  through  appropriate  modification 
of  previous  results  of  Reissner  [5]  and  Hildebrand  [6],  in  the  limiting- 
type  orthotropic  cases  for  which  Young's  modulus  in  the  direction  normal 
to  the  stressfree  edges  takes  on  the  values  zero  or  infinity,  respectively. 

THE  BOUNDARY  VALUE  PROBLEM 
We  consider  a rectangular  region  with  boundaries  x = ± a and 
y = ± c.  We  assume  that  the  boundary  portions  y = ± c are  subjected  to 

m--  ’ ' 

uniform  displacements  ±U  in  the  x-direction  and  that  displacements  in  the 
y-direction  are  prevented.  The  boundary  portions  x = ± a are  assumed 
to  be  traction  free. 

\ We  have  as  differential  equations  for  stresses  and  strains 


a + r = 

0 r + a - 0 

(1) 

* 

X,  X , y 

y.y 

1 

and 

? 

€ = 

u , 

f = V , y = U +v 

(2) 

X 

, X 

y .y  . y >x 

1 

and  therewith 

as  boundary 

conditions 

> 

y = ±c  ; 

u = ±U  , v = 0 

(3) 

x = ±a  ; 

0 = 0 , t = o 

(4) 

• 

X 

We  assume  the  medium  to  be  orthotropic  with  constitutive 

relations 

a 

a 

o a 

X 

- v-*-  . 

y - x r 

e = — — - v , y = 

(5) 

€x  ~ E 

E 

y E E y G 

X 

m 
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where  E = V EE,  with  E , E , v and  G being  constants, 
m x y x y 

For  the  case  of  isotropy  and  plane  stress  we  have 


E = E = E,  v = v and  G = 1 E/(l  + v) 

x y 


(6a) 


and  for  the  case  of  isotropy  and  plane  strain  the  corresponding  relations 


r 


. 


r i 


are 


E = E = E/(l  - vs),  v = v/(J  - v)  and  G = 1 E/(J  + v).  (6b) 

x y 


For  the  strain  energy  to  be  positive  definite  we  must  have  |v|<  1. 
The  prescribed  edge  displacements  u(x,  ±c)  = ±U  are  associated 
with  forces 

i*  a 


± P = f t(x,  ±c)dx, 
J -a 


(7) 


in  the  form 


P = KU 


(8) 


u = Uy/c,  r = U/Gc,  v = c = ct  = 0 

x y 


(9) 


with  the  corresponding  value  of  the  stiffness  coefficient  being  given  by 


K = 2Ga/c. 
o 


(9) 


UPPER  AND  LOWER  BOUND  EXPRESSIONS  FOR  K 
From  the  principles  of  minimum  potential  energy  and  minimum 
complementary  energy  we  have  as  inequalities  for  the  work  quantity  PU, 
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with  the  values  of  the  stiffness  coefficient  K being  the  principle  objectives 
of  the  following  considerations. 

As  an  elementary  approximation  we  may  assume  that  the  elastic 
region  is  in  a state  of  pure  shear,  that  is. 


1 


I = f-  fC  ra  Ie  eS  + 2vE  f f + E + G-y?l  dxdy 
d J-c  J-a  | x x mxy  y y f 


where 


a.  jf- 


•(E',  E , E ) = (E  , E , E )/( 1 - Vs ). 
x m y x m y 


For  (II)  to  be  valid  the  stress  components  a » r must  satisfy  the 
equilibrium  equations  (1)  and  the  stress  boundary  conditions  (4)  and  the 
strain  components  (,  y must  be  derived  from  displacements  uf,  "v  which 
satisfy  the  boundary  conditions  (3). 

Substituting  for  P from  (8),  equation  (1  !)  becomes  a system  of 
inequalities  for  K, 


I /U2  < K s I ,/Ua  . 
s d 


In  order  to  obtain  upper  bounds  K and  lower  bounds  K we 

u 

consider  the  set  of  displacement  and  stress  approximations 


u = U 


[c  '(c  -£)  fi({) 


£ = x/a  and  p - c/a 

and  where  primes  indicate  differentiation  with  respect  to 
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The  expressions  in  (16)  satisfy  the  displacement  boundary 
conditions  for  all  choices  of  F1  and  F2.  The  expressions  in  (20)  satisfy 
the  equilibrium  equations  (1)  for  all  choices  of  fx  and  fs . In  order  to 
satisfy  the  stress  boundary  conditions  (4)  they  must  be  such  that 

M±0  = fs(±D  = f*  (±1)  = 0 . (19) 


We  note  that  the  approximating  functions  used  in  T 1 ] are  subsets  of  the 
above  class  of  functions. 

Obvious  symmetries  pertaining  to  the  boundary  value  problem  as 
stated  require  that  Fg  and  fx  must  be  even  functions  of  £ whereas  F-  and 
f2  must  be  odd  functions  of  £. 

Using  ( 1 6) ,(  13)  and  117)  in  (12)  we  obtain,  after  completing  all 
y-integrations , 


K ,1 

u JP  / 4 , 4 

tT  ~~  2 J .{*  + T F>  + ? F 

O - 


8E*ps 

(Fa')S  +-  FgF' 


105G 


K 

K 


l 1 


vE*  4E5;c 

TE  g F=Fl  + 3Gps  Ft-  * TT <F')S} d?’ 

" 1 y \ m/ 


1 I*  1 

2 


(20) 


4 G 


3 E ps  f|  ’ 75  ~ 1571  d** 

y 


(21) 


where  = 2Ga/c.  It  remains  to  minimize  the  functional  (20)  and  to 

maximize  the  functional  (21).  We  note  that  our  elementary  approximation 

(9)  corresponds  to  the  assumption  FT  = F?  = 0 in  (16).  Furthermore, 

from  (20)  we  have  K s K . 

u o 


FIRST  APPROXIMATION  UPPER  AND  LOWER  BOUNDS 
We  obtain  a relatively  simple  formula  for  the  upper  bound  by 
setting  F3  (£)  = 0 in  the  expressions  (16)  for  u and  v.  Then  the  upper 


bound  coefficient  (20)  may  be  written  as 


IT  =2  I ,{‘  +f  r>  + 3c£  'fs  ,F>,)S}d<  • <22> 

Minimization  of  (22)  results  in  the  Euler  differential  equation 


f; 


5E* 


F,  =0 


2Gps  1 

and  the  Euler  boundary  conditions 

Fj'liD  = -f 

where,  as  stated  earlier,  Fx  (£)  is  an  odd  function  of 

The  minimum  value  of  K follows  from  (22)  and  (23)  as 


(23) 


(24) 


Kui  2 

— = 1+IF1(,) 


(?5) 


The  solution  of  (25)  subject  to  the  boundary  conditions  (24)  is, 

(26) 


sinh  Ao£/p 

F,  = - — L- 


4/i  cosh  A / P 
o o 


where 


A = J 5E*/2G 

o y 

Using  (26)  in  (25)  we  obtain  the  first  approximation  upper  bound, 


(27) 


K 

Ui 

K 


, . tanh  — 

6 Aq  P 


(28) 


We  next  calculate  a lower  bound  by  setting  (s  - 0 in  equation 
( ! 7).  The  lower  bound  coefficient  (21)  when  f?  - 0 reduces  to, 


6 


Using  (6a)  and  (6b)  wc  obtain  in  the  case  of  isotropic  plane 


stress  the  bound  inequality. 


1 - - --  P **  1 p 

^6(1  + v)  o 

and  in  the  case  of  plane  strain  the  inequality. 


(36a) 


o 


75  j\-Zv 
6 v 1 - v 


(36b) 


The  bounds  given  by  (36b)  are  identical  to  those  obtained  in  PI]. 

As  important  properties  of  the  bound  formulas  (28)  and  (34)  we  note 

that  a)  the  effect  of  the  stress  free  edges  is  to  decrease  the  stiffness 

and  b)  these  first  approximation  bound  coefficients  are  independent 

of  the  elastic  constant  E . 

x 

IMPROVED  UPPER  AND  LOWER  BOUNDS 

In  order  to  improve  the  bounds  (28)  and  (34)  we  next  consider 

the  complete  expressions  (20)  and  (21)  for  K and  K . Minimization  of 

u /, 

(20)  with  respect  to  Fx  and  F2  gives  the  Euler  differential  equations 


v 1 VEm 

Fl  2 ( ~G~ 


F'  = 0 


II  G 

2 2 E:;-'p2 

X 


f .1  -2- 

2 2 E*ps 

x 


Fi  = 0 


(37a) 


(37b) 


and  the  Euler  boundary  conditions 


f\'n>  = -f 


- vE- 

F2(1)  'I  F'*(1,=0 

x 


Appropriate  integrations  by  parts  and  use  of  (37)  and  (38)  enables 
us  to  obtain  the  minimum  value  of  (20)  in  the  form, 
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(39) 


ue  2 

iT  = 1 +T 


Upon  solving  the  system  of  differential  equations  (37)  subject  to 
the  boundary  conditions  (38)  and  substituting  for  F1  ( 1 ) in  (39)  we  obtain 


fSj.  5P  ***  - A? 

K 6 v 

o A 


Ai  A-z 

tanh  — tanh  — 
P P 


where  are  the  roots  with  positive  real  parts  of  the  equation 


- 5 F* 

- (1  .Z^,_Z  + 35G  .7 

' 10  ’ 2G  4E*  2 

L x 


m ,,  105  y 

A + = 0 

E*  4 E* 

x J x 


and  where 


f 

X 

G + vE5'* 
m 

1 

+ 

Xa(Ag  tanh 

Ap 

p 

with 

7 vE* 
m 

« 

Xi  " 6G 

Ai  As 

vE* 

m 

5E* 

I 

0 

II 

CO 

X 

2G 

Ax 

— ) + X 1 

p *3 


2 G Ai  A 2 


(1  - .7 vs)  - 


7 vE- 
m 


When  A-J  p ^ 4,  K may  be  approximated  by  its  asymptotic  form, 


ue  , J y 1 - .7  vs 

K “ P ^Aj  + A? ) / j _ 3 Ai  A2  • (44) 


Next  we  proceed  to  obtain,  in  a similar  fashion,  an  improved  lower 
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bound  for  K.  However,  maximization  of  the  functional  (21)  as  it  stands 
would  require  the  solution  of  a sixth  order  differential  equation.  In 
order  to  avoid  increasing  the  order  of  our  problem  beyond  the  fourth 
we  shall  in  what  follows,  maximize  the  expression  (21)  by  assuming 
fx  (£)  as  given  by  the  first  approximation  calculations.  In  this  way  the 
resulting  Euler  differential  equation  is  the  fourth  order  equation, 


f 


HU 

2 


42  E 525E 

Y f"  + Y f 

Gps  s E p4 
x 


1 05E 


Y 

2Gp? 


( 


m 


(45) 


with  f^  corresponding  to  fx  in  (32). 

The  constraint  conditions  on  f2  in  (19)  and  the  differential  equation 
(45)  in  conjunction  with  appropriate  integrations  by  parts  enable  us  to  write 
the  maximum  value  of  the  functional  (21)  in  the  form 


K 


is 


K 


= 1 - 


X 


tanh 


(46) 


Upon  solving  the  non-homogeneous  equation  (45)  subject  to  the 
appropriate  boundary  conditions  and  evaluating  the  integral  in  (46),  we 
obtain  as  improved  lower  bound  expression 


K 

o 


+ 


lE 

6 


X4  x(  1 + 5vG/E  ) 

o m 

>4  - (x?  + xl>x2  + xfx| 

o 1 * o 1 


(47) 


where  X.  are  the  roots  W’ith  positive  real  parts  of  the  equation 

E E 

X4  - 42  -X  + 525  = 0 , (48) 

G E 

x 

and  where 


v = x + x, 
o 1 


(49) 
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1 


with 


x = 

o 


tanh  X /p 
o 


p cosh2  X / P 


(50) 


2(X§  - X2) 

Xl  = (x*  - x|)(x*  - x?: 


(X,  tanh  x /p  - x tanh  Xi/p)(X2  tanh  X /p  - X tanh  X2/p) 

2 2 : ! 2 2 (5]) 

Xi  tanh  X2/p  - X2  tanh  X-^/p 

When  X./p  > 4,  K may  be  approximated  by  its  asymptotic  form, 

1 


K 


K 


la  _ j _£_ 

X. 


, X4  (1  + 5vG/E  )2 

i -I  ° 


m 


6 (xn  + X1)a(Xrt  + Xs): 
o 1 o 


(52) 


We  note  that  the  asymptotic  formulas  (35)  , (44),  and  (52)  correspond 
to  the  assumption  of  the  existence  of  boundary  layers  in  the  neighborhood 


of  the  edges  x = ±a.  In  the  expressions  for  K and  the  leading  terms 


represent  the  interior  solution  contribution  in  the  stiffness  coefficient 
corresponding  to  a state  of  pure  shear.  The  remaining  terms  represent 
the  boundary  layer  contributions.  Within  these  boundary  layers  the 
interior  shear  stress  U/Gc  undergoes  a rapid  transition  to  attain  the  value 
zero  at  x = ±a.  The  approximate  boundary  layer  thickness  obtained  in 
[ 1 ] for  the  isotropic  plane  strain  case  may  be  generalized  to  the  case  of 
orthotropy  by  observing  our  first  approximation  solutions.  For  the 
orthotropic  case  we  have  from  equations  (26),  (27),  (32)  and  (33)  as  a 
condition  for  the  existence  of  boundary  layers 


p 5 JE/G 


(53) 


and  the  corresponding  width  b of  these  boundary  layers  is  given  by 


b /c  = O (JG/E  ) . 
o y 


1 1 


(54) 


Our  improved  bound  calculations  show  that  there  exists  an 
additional  set  of  boundary  layers  adjacent  to  the  edges  x - ±a.  We 
demonstrate  this  by  considering  the  characteristic  equation  (48)  cor- 
responding to  the  lower  bound  calculation  which  may  be  considered  more 
relevent  in  the  context  of  a boundary  layer  estimate  pertaining  to  stress 
transitions.  As  conditions  for  the  existence  of  these  additional  boundary 
layers  we  have 

E /G  = 0(1)  , p < tffe  / E (55) 

m y x 

1 « E / G , P ^ «VE  / G ( 56 ) 

m Y 

The  widths  of  these  boundary  layers  are  given  by 

E /G  = 0(1)  , b./c  = O (t/E  /E  ) , i = 1,  2 (57) 

m i x y 

1 « E /G  , b./c  = O (JfTTg)  , b,/c  = Ol.  'G/  E ) (58) 

m 1 x **  y 

It  is  apparent  from  (54),  (57),  and  (58)  that  there  exist  two  boundary 

layers  near  each  stress-free  edge,  for  all  values  of  the  parameter  F / G. 

m 

As  far  as  the  numerical  calculation  of  the  bounds  is  concerned  we 
may  use  the  asymptotic  formulas  (35),  (44)  and  (52)  whenever, 

Min  (Re[A.,  X.].  A . X ] i 4p  . (59) 

i t o o' 


as  these  formulas  involve  the  approximation  tanh  I A,  >)/ p 1. 

The  roots  of  the  characteristic  equations  (41)  and  i 48 1 are  i _<ne  ra ! 
complex.  As  conditions  for  real  roots  we  must  have 


1 - .7  vs 
1 - v2 


+ 35 


1 - vg  Gg 

4 Eg 
m 


i v 

T 


(60) 


for  real  and  A2  and 


I for  real  and  \s  . 

The  improved  upper  (or  lower)  bound  obtained  here  is  not  valid 
when  the  equality  sign  holds  in  (60)  (or  in  (61)).  However,  expressions 
for  the  bound  coefficients  may  be  deduced  from  the  formula  (40)  (or  (47)) 
by  taking  the  limit  as  Ax  "*  As  (or  - X2). 

Isotropic  plane  stress  case.  In  the  case  of  isotropic  plane  stress  we 

use  (6a)  in  (40)  and  (47)  to  evaluate  the  improved  bounds.  Fig.  1 shows  a 

plot  of  K , K , K and  K against  p = c/a  for  the  plane  stress 
ui  us  is 

problem  when  v - 1/3.  The  plots  show  that  the  exact  formulas  and  the 
corresponding  asymptotic  formulas  (linear  functions  of  c/a)  are  in- 
^ distinguishable  in  the  range  0 < c/a  < 1,  with  c/a  = 1 representing  a 

4 square  lamina.  The  maximum  error  in  employing  the  average  of  the 

improved  bounds  in  place  of  the  exact  stiffness  coefficient  comes  out 
to  be  about  0.9%  when  c/a  = 1 and  this  error  is  linearly  decreasing 
with  decreasing  values  of  c/a.  The  corresponding  error  in  using  the 
first  approximation  bound  expressions  comes  out  to  be  about  4%. 

Orthotropic  case.  In  the  case  of  orthotropic  materials  the  bound 
s formulas  (40)  and  (47)  suggest  the  use  of  the  parameters 

^ = C/JE  E , p*  = JG/  E c/a  (62) 

x y y 

with  p representing  the  effect  of  the  elastic  coefficient  E . Fig.  2 
shows  plots  of  K , K , K and  K against  p*  when  v = 1/3,  p = 10. 

Ill  US  ^,2  J X 

As  p*  increases  the  difference  between  the  upper  bound  K and  the 

ui 

exact  result  becomes  significant,  although  the  remaining  bounds  K 

U2 

K and  K are  extremely  close.  We  note  that  for  large  values  of 
p’  t if  G and  E^  are  of  the  same  order,  c/a  » 1.  Then  our  problem  is 
equivalent  to  that  of  a vertical  beam  with  length  2c  and  thickness  2a 
subjected  to  prescribed  displacements  (zero  displacement  in  the  axial 
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and  2U  in  the  tangential  directions  at  one  end  and  the  other  end  fixed  ). 
We  have  from  (34) 


K 


■li 


2Ga 

c 


(63) 


and  when  c/a  » 1, using  the  series  expansion  of  tanh  > for  \ « 1 , we 

get 


— Ey  ag 

5 G c2 


(64) 


which  is  the  beam  stiffness  coefficient  with  a shear  deformation  factor  6/5. 


EXACT  VALUES  OF  STIFFNESS  COEFFICIENTS  FOR  TWO  LIMITING- 

TYPE  ORTHOTROPY  CASES 

Exact  expressions  for  the  stiffness  coefficients  may  be  obtained 
in  the  two  cases  of  limiting-type  orthotropy  by  following  the  formula- 
tion in  r 51  in  the  case  of  E =0  and  by  following  the  formulation  in  r61 
J x 

in  the  case  of  E = ».  These  calculations  are  presented  in  an  Appendix. 
Our  final  results  show  that 


K . = K f"l  - tanh  —1 

E = 0 o L a P-  J 


K = K 

E - oo  o 

X 


a J 

1+  P*  E 

n=0  °n 


tanh  a / p- 
n 


-i 


p *(tanh  O'  / p : ) / ry 
n n 


where 


(65) 


(66) 


a = j3/(  1 - v2)  , 


a = (n  + — ) it 
n 2 


(67) 


Fig.  3 shows  plots  of  (65)  and  (66)  for  a wide  range  of  p 
show  that  the  effect  of  E on  the  stiffness  is  small  when 


range  of  the  parameter 


x 

H £ G/JK  E . 

x y 


r’ 


These  plots 
« 1 in  the  entire 
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m/ 


m 


In  Fig.  3 we  have  also  plotted  the  mean  value  of  the  closest  bounds 

when  p = 10.  For  the  range  of  p*  shown  in  Fig.  2 the  mean  value  of 

K and  K differs  from  the  exact  values  by  less  than  3.5%  with  the 
ua  Is 

maximum  of  this  error  corresponding  to  p*  =»  4.  The  plot  in  Fig.  3 shows 
that  the  exact  result  for  E = 0 (p  = «)  may  be  used  as  an  approximation 
for  K when  p < 10,  p*  <1,  as  in  this  range  the  curves  for  p = 10  and  for 
p = 0 are  not  distinguishable. 
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APPENDIX 


a.  The  case  - 0.  We  write  the  resulting  stress  strain  relations 
in  the  form 

a = 0 , a =Ef/(l-v2),T  = Gy.  (Al) 

x y y y ' 

Proceeding  as  in  [5]  the  equilibrium  equations  (1)  are  satisfied 


by  taking 


r = f(x)  , a = -yf'(x) 

y 


The  second  equations  in  (Al)  and  (A2)  along  with  the  defining 
relation  f = v then  give  as  expression  for  v 

y .y 

1 - Vs  r 1 „ I 


- ^ ys  f ' + y + h(x) 


Satisfaction  of  the  boundary  conditions  v(x,  ±c)  = 0 gives 


g = 0 , h^c8f' 


and  therewith 


1 - v2  c2  - y2 
v = E 2 ‘ ■ 

y 


The  constitutive  relation  for  T in  (A  1 ) in  conjunction  with  (A5) 
and  the  defining  relation  y - u + v then  gives  further 

> y , x 


u _ xL  _ 1 - ya  (six  . it)  r 

" G E \ 2 6 / 


upon  taking  account  that  u is  an  odd  function  of  y. 

We  next  obtain  the  differential  equation  for  f by  setting  u(x,  c)  = U, 


1 - v8  Gc2  f = " 


I - V2 
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The  remaining  boundary  conditions  r(±  a,  y)  = 0 determine  the 
solution  of  (A7)  in  the  form. 


f=£U  2 
c 


cosh  ,/3E  / G(  1 - v2)  x/c 


cosh  ^3E  / G(  1 - v2)  a/c 

y 


Use  of  (A8)  in  the  expression  for  the  stiffness  coefficient, 
K = J*a  x('<.  c)ds/U 


gives 


'E  = 0 


G(i  - vg)  c / 3Ey  a 

3E  a an  Jg(  1 - v2)  c 

y 


b.  The  case  E^  - x . In  this  case  the  constitutive  relations  (5)  reduce 


€ -0,  f =ct/E  , y = r/G  . 

x y y y 


We  consider  u(y)  and  v(x,  y)  of  the  form 


Py  ® 
u = - i + Z u 


7±7  + I u (y)  , V = E X (x)Y  (y)  , 
ZGa  n n n 

n=0  n=0 


where  P is  a constant  to  be  determined  and,  in  order  to  satisfy  the 
conditions  v(x,  ±c)  =0, 

Yn(±c)  =0  (A  1 3) 

Equation  (A  12)  in  conjunction  with  the  constitutive  relation  for  r 
in  (All)  gives 


r = r-  + G T [ u*  + X ' Y ] , 
2a  _ n n n 

n=0 
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where  primes  and  dots  denote  differentiation  with  respect  to  the 

respective  arguments. 

From  (1)  we  have  as  an  expression  for  a > 

x 

a = - T T*  dx  . (A  1 5) 

x J y 

-a 

Use  of  (A14)  in  ( A I 5 ) and  satisfaction  of  the  boundary  conditions 

CT  (± a,  y)  = 0,  give 
x 

u”  = rx  (-a)  - X (a)lY'  /2a  . (A16) 

n n n ' n 

As  v is  odd  in  x and  even  in  y,  integration  of  (A  16)  gives 

u*  = - X (a)Y  /a  . (A  1 7) 

n n n 

From  the  second  equations  in  (All)  and  (A  1 2)  we  have 

C = E L X Y‘  . (A  1 8) 

y y n=o  n n 

Substituting  for  o from  (A18)  and  for  r from  (A14)  with  u' 
y n 

as  in  (A  17)  into  the  equilibrium  relation  between  a and  r gives 


f fE  X Y"  + GX’Y  J = 0 . 


y n n n n 


The  solutions  of  (A  19)  are  taken  in  the  form 

Y = cos  a y/c  , X = A sinh  JE  /Gey  x/c 
n n n n v y n 


where,  in  order  to  satisfy  ( A 1 3 ) , 


2n  + 1 


a = — 7" 
n 2 


and  where  the  constants  A remain  to  be  determined. 

n 
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r 


L 


We  next  use  (A  14)  to  write  for  the  boundary  values  of  r 


, » P G « I fE  [E  [E  ' 

T(a.  v)  = * 7 E v cosh  M,  — - — sinh  -Va  * 

n=0  L J G n jGnc  a \Gnc  nc 


(All) 

A formal  expansion  of  P/2a  in  a Fourier  series  using  the  ortho- 
gonal eigenfunctions  cosa  y/ c and  setting  r(a,  y)  - 0 then  gives 

n 


sinry  f !E  E /E 

n . , y a Y a . / y a 

A = P — sinh  -f-  a — - / a — cosh  i-f-  & — 

n Go  V G nc  * G nc  VG  nc 

n L 


( A23) 


From  (A14)  and  (A17),  with  Y (c)  = 0,  we  see  that 

n 


P = J t(x,  c )dx  . 


( A24) 


Integration  of  (A  17)  with  Y as  in  (A20)  and  substitution  of  the 

n 

resulting  u in  (A  12)  finally  gives 


Pc  , 2_c  > G ® 1_ 

2Ga  a J E o3 


tanh.'E  /Go  a/’c 
V n 


y n=0  n 1 -(tanh  J E /Go  a/c)/„'E  /Go  a/c 

v y n y n 


(A25) 


The  stiffness  coefficient  K„  follows  from  (A25)  as  in  (66). 

E — CD 

X 

We  note  that  the  above  calculation  is  similar  to  that  in  [6]  so  far 
as  the  solution  procedure  is  concerned.  However,  our  boundary  conditions 
are  different  from  those  considered  in  [6]. 
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Figure  1.  Values  of  stiffness  coefficient  bounds  for  isotropic  materials  with 


6h 


• 4 r 


Fig.  3.  Exact  values  of  stiffness  coefficients  Tor 


E - 0 (|i  - «),  and  R 

X X 


(p  0),  as  well  as  value 


of  (K  + K )/ 2 corresponding  to  u = F 10, 

us  12  r ' x y 

and  v 1/3  as  a function  of  p ~ JG/  E c/a. 


6 


Values  of  stiffness  coefficient  bounds  for  th 

case  p v C,/Jf.  FT  = 10  and  v = 1/3,  as  a 
x y. 

function  of  p*  a JG/V.  c/a. 
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